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This invention is best described in the following statement: 



WAVELET FILTERING 
Technical Field of the Invention 

The present invention relates generally to wavelet filters and, in particular, to 
wavelet filtering of digital images. 

Background 

The field of digital data compression and in particular digital image compression 
has attracted great interest for some time. 

In the field of digital image compression, many different techniques have been 
utilised. In particular, one popular technique is the JPEG standard, which utilises the 
discrete cosine transform to transform standard size blocks of an image into 
corresponding cosine components. In this respect, the higher frequency cosine 
components are heavily quantised so as to assist in obtaining substantial compression. 
The heavy quantisation is an example of a lossy technique of image compression. The 
JPEG standard also provides for the subsequent lossless compression of the transform 
coefficients. 

Recently, the field of wavelet transforms has gained great attention as an alternative 
form of data compression. The wavelet transform has been found to be highly suitable in 
representing data having discontinuities such as sharp edges. Such discontinuities are 
often present in image data or the like. 

Typically, data compression using wavelet techniques is a two step process. It 
comprises, firstly, a transform phase, during which the wavelet transform of the data set is 
calculated, and secondly a subsequent coding stage during which the resultant data set 
from the transform operation is separated into segments which are then coded using a 
specific coder. In decompression, the reverse occurs, with coded blocks being first 
decoded, and subsequently the inverse wavelet transform being applied to generate the 
final decompressed output. 

Typically, the wavelet transform operation utilises a number of wavelet filters. 
These filters are typically digital filters, such as finite impulse response (FIR filters) 
which ideally are linear and require a small number of taps. These filters process the 
digital images in bands and typically require a large external store for storing overlapping 
pixels, coefficients, and/or intermediate results for processing the next band. The 
bandwidth to the external store is typically limited, which restricts the rate of reading and 
writing of overlapping pixels, coefficients, and/or intermediate results to the store. 

Disclosure of the Invention 



485397 CFP1644AU Siandard09 



[I:\ELEC\CISRA\S tandard\standard09J485397.doc:pwm 




-2- 



It is an object of the present invention to substantially overcome, or at least 
ameliorate, one or more disadvantages of existing wavelet filters. 

According to a first aspect of the invention, there is provided a method of wavelet 
filtering a digital image, the digital image comprising a plurality of pixels arranged in one 
5 or more bands of pixels, the method performing the following steps for each one of said 
one or more bands of said pixels: retrieving a plurality of previous partial results from a 
local storage or remote storage; inputting a current group of adjacent said pixels; 
computing coefficients utilising said current group of adjacent pixels and said plurality of 
previous partial results; outputting said coefficients; and storing a plurality of current 
10 partial results in said local storage if said current group is not the last group in the band, in 
said remote storage if said group is the last group in the band. 

According to a second aspect of the invention, there is provided an apparatus for 
implementing the aforesaid method. 

According to a third aspect of the invention, there is provided a computer readable 
15 medium comprising a computer program for implementing the aforesaid method. 

According to a fourth aspect of the invention, there is provided a wavelet filter for 
wavelet filtering a digital image, the digital image comprising a plurality of pixels 
arranged in one or more bands of pixels, the filter comprising: a local storage for storing a 
plurality of previous partial results; a remote storage for storing a plurality of previous 
20 partial results; a controller for selecting the plurality of previous partial results from said 
local storage or said remote storage; a pixel input mechanism for inputting a current 
group of adjacent said pixels; a lifting lattice of multiplier and adder units for computing 
coefficients utilising said current group of adjacent pixels and said selected plurality of 
previous partial results; output means for outputting said coefficients; and a controller for 
25 storing a plurality of current partial results in said local storage if said current group is not 
the last group in the band, in said remote storage if said group is the last group in the 
band. 

Brief Description of the Drawings 

A number of preferred embodiments of the present invention will now be described 
30 with reference to the drawings, in which: 

Fig. 1 is a schematic representation of a K-line 5-3 wavelet filter (with K = 8) for 
filtering a digital image in accordance with a first embodiment; 

Fig. 2 is a schematic representation of the cascading of a wavelet filter of Fig. 1 in 
time; and 
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Fig. 3 is a flow diagram of a method for use in wavelet filtering in accordance with 
a second embodiment; 

Fig. 4 is a schematic representation of a general-purpose computer for use in 
implementing the method of Fig. 3. 
5 Detailed Description including Best Mode 

Where reference is made in any one or more of the accompanying drawings to steps 
and/or features, which have the same reference numerals, those steps and/or features have 
for the purposes of this description the same function(s) or operation(s), unless the 
contrary intention appears. 

10 The principles of the preferred embodiments described herein have general 

applicability to the implementation of K-line wavelet filters. However, for ease of 
explanation, the embodiments has been described with reference to a K-line 5-3 wavelet 
filter for use in a forward discrete wavelet transform of digital images. However, it will be 
readily evident that the invention is not limited thereto. For instance, the wavelet filter can 

15 also be used in K-line 9-7 wavelet filters. Furthermore, it would be apparent to a person 
skilled in the art, that the wavelet filter can be modified for use in an inverse discrete 
wavelet transform. Still further, it would be apparent to a person skilled in the art that two 
such wavelet filters can be combined in known manner to provide a 2-dimensional 
discrete wavelet transform or modified for use in a 2-dimensional inverse discrete wavelet 

20 transform. Moreover, the wavelet filter can be used in many different applications. For 
examples of the many different applications of wavelet analysis to signals, reference is 
made to a survey article entitled "Wavelet Analysis" by Bruce ET. AL appearing in IEEE 
spectrum, October 1996 pages 26 to 25. For a discussion of the different applications of 
wavelets in computer graphics, reference is made to "Wavelets for Computer Graphics", 

25 page 5, 1. Stollinitz et. al. published 1996 by Morgan Kaufnann Publishers, Inc. 

As mentioned above, the embodiments of the invention can be used in a discrete 
wavelet filter transform as well as an inverse discrete wavelet transform. For sake of 
simplicity, the term wavelet transform used herein or variations thereof, is taken to 
include both a forward and/or an inverse wavelet transform, unless the contrary intention 

30 appears. Similarly, the term pixel as used herein or variations thereof, is taken to taken to 
refer to the original pixels or wavelet coefficients of a digital image, unless the contrary 
intention appears. Similarly, the term digital image or variations thereof is taken include 
an original image or a sub-band of that original image, or one or more associated sub- 
bands of wavelet transformed coefficients of that original image. 



485397 CFP1644AU Standard09 



[I:\ELEC\CISRA\S tandard\standard09]485397.doc:pwm 




-4- 



Turning now to Fig. 1, there is shown a schematic representation of a K-line 5-3 
wavelet filter (with K = 8) for filtering a digital image in one dimension in accordance 
with a first embodiment. The digital image comprises a plurality of pixels arranged in a 
plurality of columns in the vertical direction and a plurality of rows in the horizontal 
5 direction. 

The filter 100 processes one horizontal band of the digital image after another. The 
digital image comprises one or more such horizontal bands, with the horizontal bands 
comprising a plurality of said rows of said pixels. The filter 100 initially commences 
processing the uppermost band of the image, and then processes the next adjacent band of 
10 the image, and so on until the lowermost band of the image is processed. Within each 
band, the filter 200 processes one column of pixels after another. The filter 100 first 
processes the left column of the band, and then the next adjacent column of the band, and 
so on until the right column of the band is reached. For ease of understanding, the pixels 
in each column of a band from the topmost pixel to the lowermost pixel are represented as 

15 XQ t Xj t X2, . . . XK'-l- 

Within each column of a band, the filter 100 processes a first group of adjacent 
pixels (commencing with the topmost group) and then processes a second group of 
adjacent pixels adjacent the first group. 

In this way, the filter 100 performs a one dimensional discrete wavelet transform in 

20 the vertical direction. In another embodiment, the filter 100 can also be adapted to 
perform a one dimensional discrete wavelet transform in the horizontal direction by 
processing vertical bands one after another in similar fashion. In still another 
embodiment, the filter can also be adapted to provide a 2-dimensional wavelet transform. 
In the latter case, one filter 100 processes one column at a time first, and a second filter 

25 100 after that collects the output of the first filter 100 and processes one row at a time. 

The pixels ( e.g. xo t x/ t x 2i x 3i x 4t xs t x^, x 7t xg t xg t x/o, xn, x/j. x/4 and x/5) of a column 
of a band are fed to the filter 100. These pixels are fed in parallel to the wavelet filter 100 
via a number 2K (with K = 8) of input lines 101. The 16 input lines 101 are connected to 
eight 2-input multiplexers 103-1, 103-2, 103-3, 103-4, 103-5, 103-6, 103-7, and 103-8. 

30 The multiplexers can select either the bottom eight lines or top eight lines for input of the 
groups of pixels (x 0t x lt x 2i x 3t x 4t x 5> x 6t and x 7 ) or {x 8t x 9t x Wt x 1!t x !2 , x /3i x /4 and x I5 ) 
respectively. The multiplexers 103-1, 103-2, 103-3, 103-4, 103-5, 103-6, 103-7, and 103- 
8 are connected to a lifting lattice of multiplier and adder units 102, 104, 106, 108, 110, 
112, 114, and 116, which compute the high pass coefficients and low pass coefficients on 
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the top 8 lines (ie. coefficients (x$. jcp. x/o. xu_ xn, x/3. x/ 4 and x/5)) or the bottom 8 lines ( ie. 
coefficients (xo. *?. X4, X5. X6. and X7)). The lifting lattice of multiplier and adder units 
102, 104, 106, 108, 1 10, 1 12, 1 14, and 1 16, perform the following calculations: 

5 d n =x 2n +a{x 2n _ 1 +x 2n+i ) n = 0,l,...,K-l (1) 

and, 

s„ = *2„ +I + P{d n + d n+1 ) n = -1,0,..* - 2 (2) 

where d n are all the high pass coefficients, s n all the low pass coefficients, a = -0.5, and P 
10 =0.25. 

For example, the multiplier and adder units 102, 104 106, 108, 1 10, 1 12, 1 14, and 
1 16 perform the following calculations on the selected input pixels (xq, *2, *3, *4, *5 t *6, 
and X7) to compute the high pass coefficients d n and the low pass coefficients s n \ 



unit 102 


d 0 =x 0 +a(x_ l + x x ) 


(3) 


unit 104 


d x =x 2 + a{x x + x 3 ) 


(4) 


unit 106 


d 2 = x 4 + a(x 3 + x 5 ) 


(5) 


unit 108 


d^ = x 6 + a{x 5 + x 7 ) 


(6) 


unit 110 


= + + d-0 ) 


(7) 


unit 112 




(8) 


unit 114 


s x = x 3 +/?(</, + d 2 ) 


(9) 


unit 116 


S 2 = X 5 4- /?(^ 2 + ^ 3 ) 


(10) 



The high pass and low pass coefficients d n and s n are then passed to other discrete 
25 wavelet transform circuits or entropy coder (not shown in the figure). 

The multiplier and adder units 1 10 and 102 also accept inputs from either the partial 
result registers 109 or external buffer 107, depending on the control on the multiplexer 
105. Namely, the multiplier and adder units 110 and 102 accept input from the external 
buffer 107 when computing the first 8 lines and accepts input from the partial result 
30 registers 109 for the last 8 lines. 

A controller 1 1 1 controls the storing of a last input pixel and a last high pass 
coefficient (e.g. x 7 and d$ ) to the partial result register 109 or external buffer 107 
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depending on which lines are being processed. Namely, the last input pixel and the last 
high pass coefficient of the first 8 lines are stored in the partial result register 109. 
However, when the lifting lattice of multiplier and adder units is computing coefficients 
for the last 8 lines, the last input pixel and last high pass coefficient are stored in the 
5 external buffer 107. The last input pixel and last high pass coefficient are required for 
computation of the low and high pass coefficients of the next adjacent group of pixels 
within the column. 

For example, in the case where the filter 100 is currently processing the pixels (x 8 . 
x 9t xjo. xu, x/2, xj3, xj4 , and x/ 5 ), the multiplier and adder units 110 and 102 accept as input 

10 the last pixel and coefficient x 7 and d 3 from the partial result registers 109, which were 
previously stored in the partial result register 109 during the previous computation of the 
pixels (x 0t x it x 2 , x 3i x 4t x 5i X6, and x 7 ) of the same column and band. In the case where the 
filter 100 is currently processing the pixels (x 0 to x 7 ), the multiplier and adder units 110 
and 102 accept as input the last pixel and coefficient x I5 and d 7 from the external buffer 

15 107. These inputs x }5 and d 7 were previously stored in the external buffer 107 during the 
previous computation of the pixels (x 8 to x 15 ) of the same column but in the last previous 
band. In the case where there are no previous pixels (e.g. topmost band) the boundary is 
symmetrically reflected. 

Figure 2 illustrates the cascading of the wavelet filter of Fig. 1 in time. In cycle 0, 

20 for the first set of inputs the wavelet filter 100 either symmetrically reflects the pixels 
near the boundary or reads the required partial results from external buffer 107. The filter 
100 then stores the partial results into the partial results register 109, and outputs the low 
and high pass coefficients. In cycle 1, the lifting lattice is 'moved' to the next set of input 
lines, reads input from partial results register 109, and outputs the next low and high pass 

25 coefficients. The partial results are stored in the partial results register 109 and external 
buffer 107. 

The first embodiment has been described with respect to a filter 100 having two 
groups of eight input lines 101 and eight two-input multiplexers 103-1 to 103-8. It would 
be apparent to a person skilled in the art, that the number of groups can be increased (e.g. 
30 N = 4 groups) with a corresponding modification of the eight multiplexers, (ie. eight N- 
input multiplexers). The external buffer 107 is a slow memory storage and requires at 
* least a memory size equivalent to the size of an input pixel and a high pass coefficient 
times the number of columns. The filter 100 stores in external memory the last input pixel 
and the last high pass coefficient x /5 and d 7 of each column of each band for further 
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processing in the next band. On the other hand, the partial results register is a fast local 
storage and requires only a memory size of an input pixel and high pass coefficient. 

The effective value of K' (the width of the band) can be increased by increasing the 
number of groups of pixels (eg. two groups, K' = 2xK = 16) or (four groups K' = 4xK 
5 =32)). In the present embodiments, this has the advantage of reducing the number of 
accesses to the slow external memory. Namely, where there are multiple groups, the last 
pixel and last high pass coefficient can be stored in the fast local memory, namely the 
partial results register 109, rather than the slow external memory 107, thus enabling the 
wavelet transform to be performed more efficiently. Of course, the last pixel and last high 

10 pass coefficient of the last group is required to be stored in external memory. 
Furthermore, the amount of hardware required to implement the filter is not significantly 
increased. Still further, the filter has the capability of filtering with different values of K\ 
In an alternate embodiment, the multiplexers 103-1, 103-2, 103-3, 103-4, 103-5, 
103-6, 103-7, and 103-8 may be dispensed with. In this case, the multiplexing is done by 

15 pixel fetching mechanism implicitly. For instance, the pixels are fed to the lifting lattice 
from a frame buffer by an address generator. 

Turning now to Fig. 3, there is shown a method for use in wavelet filtering in one 
dimension in accordance with a second embodiment. The method comprises in part a sub- 
procedure 300 for use in the wavelet filtering of a digital image. As mentioned 

20 previously, the digital image comprises a plurality of pixels arranged in a plurality of 
columns in the vertical direction and a plurality of rows in the horizontal direction. The 
digital image comprises one or more horizontal bands, each of which comprise a plurality 
of said rows of pixels. The method processes one horizontal band of the digital image 
after another. The method initially commences processing the uppermost band of the 

25 image, and then processes the next adjacent band of the image, and so on until the 
lowermost band of the image is processed. Within each band, the method processes one 
column of pixels after another. The method first processes the left column of the band, 
and then the next adjacent column of the band, and so on until the right column of the 
band is reached. For ease of understanding, the pixels in each column of a band from the 

30 topmost pixel to the lowermost pixel are represented as x 0i x It x 2i ... x K >. h 

The method of wavelet filtering calls the sub-procedure 300 for processing each 
column of a band. Within each column of a band, the sub-procedure 300 processes in turn 
a first group of K adjacent pixels, a second group of K adjacent pixels adjacent the first 
group , a third group of K adjacent pixels adjacent the second group, and so on to the last 
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group of pixels. The top and bottom pixel of the first and last group of pixels respectively 
being located adjacent the boundary of the band. 

In this way, the method of wavelet filtering performs a one dimensional discrete 
wavelet transform in the vertical direction. In another embodiment, the method can also 
5 be adapted to perform a one dimensional discrete wavelet transform in the horizontal 
direction by processing vertical bands one after another in similar fashion. In still another 
embodiment, the method can also be adapted to provide a 2-dimensional wavelet 
transform. In the latter case, the method processes one column at a time first, and the 
method after that collects the output and processes one row at a time. 

10 The sub-procedure commences at step 301 and proceeds to step 302 where a 

counter M is set to one. As mentioned previously, the sub-procedure processes the pixels 
within a column of a band. The sub-procedure continues to a loop 306, 308, 310, 312, 
314, and 316, where groups of a number K adjacent pixels are retrieved and processed 
one after another. A first group of K adjacent pixels (nearest the top boundary) is 

15 processed during the first pass of the loop 304 to 313, a next group of K adjacent pixels 
(adjacent the first group) is processed during the second pass of the loop 304 to 313, and 
so on. The processing of one group of pixels after another is achieved by incrementing the 
counter M. The sub-procedure terminates once K' pixels have been processed, when M.K 
= K\ Preferably, K = 4 and K'= 4 or 8. 

20 After the counter M has been set to one, the sub-procedure proceeds to step 304. In 

step 304, the method reads the corresponding previous partial results from an external 
buffer. These previous partial results having being stored in the external buffer during the 
processing of the sub-procedure of the same column in the last previous band of pixels. In 
the case where there are no previous pixels (e.g. topmost band) the boundary is 

25 symmetrically reflected. 

After step 304, the method proceeds to loop 304 to 316. This loop will be 
described with respect to the general case where M = M\ During step 306, the sub- 
procedure reads the M'th group of K adjacent pixels and then computes the low and high 
pass coefficients in step 308 in accordance with the equations (1) and (2) described above. 

30 In computing these high and low pass coefficients the method takes as input the M'th 
group of K adjacent pixels. It also accepts input from either the external buffer or a local 
register which stores the last pixel and the last high pass coefficient. The step 308 outputs 
to other discrete wavelet transform or entropy encoding procedures (not shown in the 
figure). 
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For example, in the case where the sub-procedure is currently processing the pixels 
(x s , x 9t x/o, x !it xj2, xi 3t x/4 , and x I5 ) 9 the sub-procedure step 308 accepts as input the last 
pixel and high pass coefficient x 7 and d 3 from the partial result registers, which were 
previously stored in the partial result register during the previous computation of the 
5 pixels (x 0i x it X2 t X3, x 4 , x 5i x 6t and X7) of the same column and band. In the case where the 
sub-procedure is currently processing the pixels (xo to x 7 ) y the sub-procedure step 308 
accepts as input the last pixel and high pass coefficient from the external buffer. These 
latter inputs were previously stored in the external buffer during the previous computation 
of the last group of pixels of the same column but in the last previous band. In the case 
10 where there are no previous pixels (e.g. topmost band) the boundary is symmetrically 
reflected. 

Ater step 308, the method proceeds to a decision block 310, where a check is made 
whether M.K = K\ If the decision block returns false, then the method proceeds to step 
312, where the last pixel and last high pass coefficient are stored in the local partial 
15 results register for use in the computation 308 of the next group of pixels. 

In the next step 314, the counter M is incremented and the last pixel and the last 
high pass coefficient are read 316. The method then returns to step 306 for the next pass 
of the loop 306 to 3 16. In the event, the decision block 310 returns true, then the last pixel 
and the last high pass coefficient are stored in the external buffer and the sub-procedure 
20 finishes for the current column of the band. 

As can be seen, the loop can be repeated as many times as desired to increase the 
value of K\ 

The aforementioned preferred method comprise a particular control flow. There are 
many other variants of the preferred method which use different control flows without 

25 departing the spirit or scope of the invention. Furthermore one or more of the steps of the 
preferred method may be performed in parallel rather sequential. 

The method of Fig. 3 is preferably practiced in dedicated hardware, such as shown 
in Fig. 1. However, the method of Fig. 3 can also be practiced using a conventional 
general-purpose computer system 400, such as that shown in Fig. 4 wherein the processes 

30 of Fig. 3 may be implemented as software code, such as an application program executing 
within the computer system 400. In particular, the steps of method of Fig. 3 are effected 
by code instructions in the software that are carried out by the computer. This method has 
particular application in computer systems restricted to large slow external buffers and 
small fast local memories. 
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The software may be stored in a computer readable medium, including the storage 
devices described below, for example. The software is loaded into the computer from the 
computer readable medium, and then executed by the computer 

The use of the computer program readable medium together with the software in 
the computer system preferably effects an advantageous apparatus in accordance with the 
embodiments of the invention. 

The computer system 400 comprises a computer module 401, input devices such as 
a keyboard 402 and mouse 403, output devices including a printer 415 and a display 
device 414. A Modulator-Demodulator (Modem) transceiver device 416 is used by the 
computer module 401 for communicating to and from a communications network 420, for 
example connectable via a telephone line 421 or other functional medium. The 
modem 416 can be used to obtain access to the Internet, and other network systems, such 
as a Local Area Network (LAN) or a Wide Area Network (WAN). 

The computer module 401 typically includes at least one processor unit 405, a 
memory unit 406, for example formed from semiconductor random access memory 
(RAM) and read only memory (ROM), input/output (I/O) interfaces including a video 
interface 407, and an I/O interface 413 for the keyboard 402 and mouse 403 and 
optionally a joystick (not illustrated), and an interface 408 for the modem 416. A storage 
device 409 is provided and typically includes a hard disk drive 410 and a floppy disk 
drive 411. A magnetic tape drive (not illustrated) may also be used. A CD-ROM 
drive 412 is typically provided as a non- volatile source of data. The components 405 
to 413 of the computer module 401, typically communicate via an interconnected bus 404 
and in a manner which results in a conventional mode of operation of the computer 
system 400 known to those in the relevant art. Examples of computers on which the 
embodiments can be practised include IBM-PC's and compatibles, Sun Sparcstations or 
alike computer systems evolved therefrom. 

Typically, the application program of the preferred embodiment is resident on the 
hard disk drive 410 and read and controlled in its execution by the processor 405. 
Intermediate storage of the program and any data fetched from the network 420 may be 
accomplished using the semiconductor memory 406, possibly in concert with the hard 
disk drive 410. In some instances, the application program may be supplied to the user 
encoded on a CD-ROM or floppy disk and read via the corresponding drive 412 or 411, 
or alternatively may be read by the user from the network 420 via the modem device 416. 
Still further, the software can also be loaded into the computer system 400 from other 
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computer readable medium including magnetic tape, a ROM or integrated circuit, a 
magneto-optical disk, a radio or infra-red transmission channel between the computer 
module 401 and another device, a computer readable card such as a PCMCIA card, and 
the Internet and Intranets including email transmissions and information recorded on 
websites and the like. The foregoing is merely exemplary of relevant computer readable 
mediums. Other computer readable mediums may be practiced without departing from 
the scope and spirit of the invention. 

The aforementioned embodiment shown in Fig. 1 can be modified for use in a 5/3 
inverse transform wavelet filter. In this modified embodiment, the 5/3 inverse transform 
wavelet filter calculates the original pixels from the low pass coefficients and high pass 
coefficients in accordance with the following formulae: 

P(d n _,+d n ) n = 0.1 K-l (11) 

and, 

X 2n+1 = d n-a(x2n+X2n + 2) * = "UU^ " 2 (12) 

where a = -0.5 and p = 0.25. The multiplier and adder units of Fig. 1 will be modified 
accordingly to compute the pixels in accordance with equations (11) and (12). This filter 
takes as input the coefficients s 0 to Sk-i and d 0 to d k -i , and also requires as input the 
coefficients d.j and x. 2 from the previous band, which are stored in the partial register 109 
or external buffer 107. 

The aforementioned embodiment of Fig. 1 can also be modified without departing 
from the spirit of the invention to other wavelet filters such as a 9/7 wavelet filters. 

In this modified embodiment, the 9/7 forward transform wavelet filter computes the 
low pass coefficients and high pass coefficients in accordance with the following 
formulae: 

< I =* 2 n + t+ a(x 2n +x 2n+2 ) n = 0,1,...,K-1 (14) 

Sn = *2n + P (d „_ 7 + d ' n ) tl = 0,1 K - 1 (15) 

d n =d' n +y(s' n + s n+1 ) n= - 1,0 K - 2 (16) 

s n =s; ^Sid^ +dj n = -1,0,...,K- 2 (17) 
where d n are all the high pass coefficients, s n all the low pass coefficients, d' n and s '„ are 
intermediate values, a = -1.5861, (3 = -0.052980, y = 0.88291, and 5= 0.44351. The 
multiplier and adder units of Fig. 1 will be modified accordingly to compute the 
coefficients in accordance with equations (14) to (17). This filter takes as input pixels x f 
to Xk and also requires as input the pixel x 0 , the intermediate values d'./ and s'./, and the 
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high pass coefficient d. 2 from the previous band, which are stored in the partial register 

109 or external buffer 107. 

In another modified embodiment, the 9/7 inverse transform wavelet filter computes 

the original pixels from the low pass coefficients and high pass coefficients in accordance 

with the following formulae: 

s n =s n ~ S(d n ^d n ) n = 0,1,...,K-1 (18) 
d' n =d n - y(s n + s n+1 ) n = - 1,0,...,K - 2 (19) 

x 2n =s n ~fi(d n ^d n ) n=-1,0,..-,K-2 (20) 

X*i + 1 = d 'n -<*( * 2 n +*2n + 2 ) " = "2,-f , K - 3 (21) 

where d n are all the high pass coefficients, s n all the low pass coefficients, d ' n and s * n are 
intermediate values, a = -1.5861, p = -0.052980, y - 0.88291, and 8= 0.44351. The 
multiplier and adder units of fig. 1 will be modified accordingly to compute the pixels in 
accordance with equations (18) to (21). This filter takes as input coefficients sqXo Sk-i and 
d 0 to d k -i and also requires as input the pixel x. 4 , the intermediate values d\ 2 and s '.y, and 
the high pass coefficient d.j from the previous band which are stored in the partial 
register 109 or external buffer 107. 

The aforementioned embodiment shown in Fig. 3 can also be modified in similar 
fashion to other wavelet filters, such as a 5/3 inverse transform wavelet filter or 9/7 
wavelet filters. 

Industrial Applicability 

It is apparent from the above that the embodiment(s) of the invention are applicable 
to the computer and data processing industries. In particular, the embodiments have 
application in digital data compression of digital images and videos. 

The foregoing describes only one embodiment/some embodiments of the present 
invention, and modifications and/or changes can be made thereto without departing from 
the scope and spirit of the invention, the embodiment(s) being illustrative and not 
restrictive. 

In the context of this specification, the word "comprising" means "including 
principally but not necessarily solely" or "having" or "including" and not "consisting only 
of. Variations of the word comprising, such as "comprise" and "comprises" have 
corresponding meanings. 
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The claims defining the invention are as follows: 

1 . A method of wavelet filtering a digital image, the digital image comprising a 
plurality of pixels arranged in one or more bands of pixels, the method performing the 

5 following steps for each one of said one or more bands of said pixels: 

retrieving a plurality of previous partial results from a local storage or remote 
storage; 

inputting a current group of adjacent said pixels; 

computing coefficients utilising said current group of adjacent pixels and said 
10 plurality of previous partial results; 

outputting said coefficients; and 

storing a plurality of current partial results in said local storage if said current group 
is not the last group in the band, in said remote storage if said group is the last group in 
the band. 

15 

2. A method as claimed in claim 1, wherein said digital image is an original image and 
said wavelet filtering performs a forward wavelet transform. 

3. A method as claimed in claim 1, wherein said digital image is a sub-band of an 
20 original image and said pixels are coefficients and said wavelet filtering performs a 

forward wavelet transform. 

4. A method as claimed in claim 1, wherein said digital image is one or more 
. associated sub-bands of an original image, and said pixels are coefficients and said 

25 wavelet filtering performs an inverse wavelet transform. 

5. A method as claimed in claim 1, wherein said method further comprises: 
repeating said retrieving, inputting, computing, outputting, and storing steps in 

sequence a plurality of times. 

30 

6. A method as claimed in claim 2 or 3, wherein said computing step is calculated in 
accordance with a 5/3 forward wavelet transform. 
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7. A method as claimed in claim 4, wherein said computing step is calculated in 
accordance with a 5/3 inverse wavelet transform. 

8. A method as claimed in claim 2 or 3, wherein said computing step is calculated in 
5 accordance with a 9/7 forward wavelet transform. 

9. A method as claimed in claim 4, wherein said computing step is calculated in 
accordance with a 9/7 inverse wavelet transform. 

10 10. A method as claimed in claim 1, wherein said plurality of partial results stored in 
said local storage comprises a pixel and a coefficient. 

11. A method as claimed in claim 1 , wherein said plurality of partial results stored in 
said local storage comprises a pixel, intermediate values, and a coefficient. 

15 

12. Apparatus for implementing the method as claimed in any one of the preceding 
claims. 

13. A computer readable medium comprising a computer program for implementing the 
20 method as claimed in any one of the preceding claims. 

14. A wavelet filter for wavelet filtering a digital image, the digital image comprising a 
plurality of pixels arranged in one or more bands of pixels, the filter comprising: 

a local storage for storing a plurality of previous partial results; 
25 a remote storage for storing a plurality of previous partial results; 

a controller for selecting the plurality of previous partial results from said local 
storage or said remote storage; 

a pixel input mechanism for inputting a current group of adjacent said pixels; 

a lifting lattice of multiplier and adder units for computing coefficients utilising said 
30 current group of adjacent pixels and said selected plurality of previous partial results; 

output means for outputting said coefficients; and 

a controller for storing a plurality of current partial results in said local storage if 
said current group is not the last group in the band, in said remote storage if said group is 
the last group in the band. 
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15. A wavelet filter as claimed in claim 12, wherein said pixel input means comprises a 
plurality of multiplexers for selecting in turn a group of pixels as said current group. 

16. A wavelet filter as claimed in claim 12, wherein said digital image is an original 
image and said wavelet filtering performs a forward wavelet transform. 

17. A wavelet filter as claimed in claim 12, wherein said digital image is a sub-band of 
an original image and said pixels are coefficients and said wavelet filtering performs a 
forward wavelet transform. 

18. A wavelet filter as claimed in claim 12, wherein said digital image is one or more 
associated sub-bands of an original image, and said pixels are coefficients and said 
wavelet filtering performs an inverse wavelet. 

19. A wavelet filter substantially as described herein with reference to the 
accompanying drawings. 

20. A method of wavelet filtering, the method substantially as described herein with 
reference to the accompanying drawings. 

21. A computer readable medium comprising a computer program, the computer 
program substantially as described herein with reference to the accompanying drawings. 
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Abstract 



WAVELET FILTERING 

The preferred method wavelet filters a digital image. The digital image comprises a 
plurality of pixels arranged in a plurality of columns and rows. The method processes one 
or more bands of pixels of the image in turn, where each band comprises a plurality of 
rows. The method performs the following steps for each column in each one of the one or 
more bands of pixels. The method retrieves (316) a last pixel and last high pass 
coefficient previously stored in local storage. The method also inputs (306) a current 
group of adjacent pixels. The method then computes (308) the low and high pass wavelet 
coefficients utlising the current group of adjacent pixels and the previous last pixel and 
last high pass coefficient. The method then outputs the low and high pass coefficients, 
and stores (312, 318) a current last pixel and high pass coefficient in either the local 
storage or external memory. 

Fig. 3 
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